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Abstract 

We construct SU(n + 1) BPS spherically symmetric monopoles with minimal sym- 
metry breaking by solving the full Weyl equation. In this context, we explore and 
discuss the existence of open spin chain-like part within the Weyl equation. For in- 
stance, in the SU(3) case the relevant spin chain is the 2-site spin 1/2 XXX chain 
with open boundary conditions. We exploit the existence of such a spin chain part in 
order to solve the full Weyl equation. 



1 Introduction 



The integrability of the self-dual Yang-Mills equation was first realized by Ward [2], who 
demonstrated that the twistor transform of Penrose could be used to provide a correspon- 
dence between instantons and certain holomorphic vector bundles over the twistor space CP 3 
[3] . Since then many alternative methods (see jl] and References therein) have been proposed 
in order to find instanton and monopole solutions of the self-dual Yang-Mills equation. 

A powerful approach introduced by Nahm [5[ [6j is the so-called Atiyah-Drinfeld-Hitchin- 
Manin-Nahm (ADHMN) construction. The Atiyah-Drinfeld-Hitchin-Manin (ADHM) ap- 
proach [7] , allows the construction of instantons in terms of linear algebras in a vector space, 
which dimension is related to the instanton number. Since monopoles correspond to infinite 
action instantons, an adaption of the ADHM construction involving an infinite dimensional 
vector space might also be possible. Nahm was able to formulate such an adaption in the 
ADHMN construction. To perform this construction, a nonlinear ordinary differential equa- 
tion (i.e. the Nahm equation) must be solved and its solutions (i.e. the Nahm data) are 
used to define the Weyl equation. 

The Nahm equations provide a system of non-linear ordinary differential equations for 
three nxn anti-hermitian matrices Tj (the so-called Nahm data) functions of the variable s: 



where n is the magnetic charge of the BPS monopole configuration and e^k is the totally 
antisymmetric tensor. 

In the ADHMN approach, the construction of SU(n + 1) monopole solutions of the Bo- 
gomolny equation with topological charge n is translated to the following problem, known as 
the inverse Nahm transform [5]. Given the Nahm data for a n-monopole the one-dimensional 
Weyl equation 



for the complex In- vector v(x, s), must be solved. I„ denotes the nxn identity matrix, 
x = (xi,x 2 ,x^) is the position in space at which the monopole fields are to be calculated 
and (<7i, cr 2 , (T3) are the familiar 2x2 Pauli matrices. In the minimal symmetry breaking, 
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dTi _ 1 
~~ds ~ 2 
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the Nahm data T/s can be cast as 



Ti 2 fi 



I = 



1,2,3 



(3) 



where Tj's form the n- dimensional representation of su(2) and satisfy 



(4) 



Let us choose an orthonormal basis for these solutions, satisfying 



/ 



v' v ds = I. 



(5) 



Given v(x, s), the normalized vector computed from (|2J) satisfying dSJ), the Higgs field $ and 
gauge potential A k are given by 



where the integrations are to be performed over the range spanned by the minimum and 
maximum eigenvalues of the asymptotic form of $. Then the corresponding Higgs field and 
gauge potentials satisfy the self-dual Yang-Mills equations, and they are smooth functions 
of x. 

Recently in [1J, SU(n+l) (for generic values of n) spherically symmetric BPS monopoles 
with minimal symmetry breaking located at the origin were obtained by solving the Weyl 
equation for x = (0,0, r) (see also [8] for generic spherical monopoles). In this case, the 
corresponding fields Q were calculated along the axis (0,0, r). However, their azimuthal 
dependence can not be implemented in a straightforward manner. To do so either the full 
Weyl equation for x = (x, y, z) should be solved or, a suitable similarity transformation 
should be introduced in order to derive the corresponding solutions in terms of the simple 
ones of pp. That way all the gauge fields can be recovered. 

As will be transparent in the next section, solving the full Weyl equation for the generic 
case is a rather technically difficult task. However, the crucial observation is the existence of 
a spin chain-like part within the full Weyl equation. This observation allows us to bypass the 




(6) 
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existing technical difficulties. More precisely, we present a new approach for constructing the 
spherically symmetric monopole solutions of the self-dual Yang-Mills equation by connecting 
the corresponding full Weyl equation with a one- dimensional quantum integrable model, that 
is the XXX spin chain (also knows as, isotropic Heisenberg model). To our knowledge this is 
the first attempt to make such a connection; although, the integrability of the Weyl equation 
was also discussed in [9]. 

In fact, it is worth noting that the effective Yang-Mills dynamics in several special lim- 
its is described by completely integrable systems related to Heisenberg spin chain and its 
generalization (see, e.g. [10] and References within). It was shown for instance in that 
the one loop mixing matrix for anomalous dimensions in N = 4 super-symmetric Yang-Mills 
theory can be identified with the conventional isotropic Heisenberg spin 1/2 chain. Also, it 
was shown in [T2"l [13] that the XXX Heisenberg chain of noncompact spin s — 0, describes 
the high energy scattering of hadrons in multi-color Quantum Chromo-Dynamics [12] . 

2 BPS Monopoles 

2.1 577(2) BPS Monopoles 

Recall the full Weyl equation (T5]) takes the form 



d_ 

ds 



+ iTj <g> Oj — Xj <g> Oj ) v(x, s) = 0, (7) 



where (x\ } X2, £3) = (x, y, z) and the vector v(x, s) is of the form v 

The spin zero representation of the su(2) case is trivial (that is, Tj = 0). Therefore, the 
corresponding full Weyl equation ([7]) for the SU{2) BPS monopole reduces to: 





where we define u = cos 6, ip = sin^e and its complex conjugate ip and thus, 

z = ru, x — iy = rip, x + iy = rip. (9) 
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Notice that u 2 + \ip\ 2 = 1. Then, equation (JSJ) reduces to the following set of differential 
equations 

Vi — ru v\ — rip Vi — 

ii 2 + ru v 2 — rip v 1 = 0, (10) 

and its solution is given by 

Vi = K\ e rs + k 2 e~ rs 

V2 = ^(l-u)e rs -^(l + u)e- rs . (11) 

Next choose a vector v to be of the form (TTTT) with constants R\, k 2 , such that the set {v, v} 
defines an orthonormal system. Thus, introducing the scalar product 

<v,v>= J v\s)v(s)ds, (12) 

the following relations must be satisfied 

< v, v > = < v, v >= 1, <v,v> = 0. (13) 

Then, due to the constraints ([TBI , the solutions ffTTj) lead to the following choice of constants 



K\ — \l (1 + u), K 2 = \ (1 — u) 

l/ 4sinh2r v ; ' V 4sinh2r v ; 



" 1 = V4lmh2^^ " 2 = -Viimh2^^ (14) 

Therefore, we end up with two orthonormal vectors of the form 



4 sinh 2r 



[1 + u) e rs + (1 -u) e~ rs 
■0 e rs — ip e~ rs 

ip e rs — ip e~ rs 
1-u) e rs + (l + u) e~ rs 



(15) 



4 sinh 2r 

Any other choice of constants is gauge equivalent to the above ones (due to gauge freedom). 
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Let v = vi and v = V2. Then, the associated Higgs field and gauge potentials are 2x2 
matrices with elements 

$jj = — i J s vjvj ds 

A k .. = J vj d k Vj ds. (16) 



2.2 SU(n + 1) BPS Monopoles 



The n- dimensional representation of su(2) is of the form 

n— 1 n— 1 n 



n = X) (4l+i + e Sifc) ' r2 = i E ft ( e Sit ~ e **+i) > T 3 = X ° fc e *fc } ( 17 ) 

fe=l A;=l fc=l 

where e ^ are n x n matrices defined by: ( ) = 8i k Sji and 

1 \ / kl 

a k = n + l-2k, C k = y/k(n-k). (18) 

The Nahm data for the SU(n+l) spherically symmetric monopoles of charge n are given 
by ([3]) where fi — f — —-■ Assume that the vector v(x, s) is of the form 

n 

v(x, s) = hi(r, s) t } ® (<7i(r, s) eg + <? 2 (r, s) eg) (19) 
1=1 

where ej^ is the n-dimensional column vector with one at the position k G Z + and zero 
elsewhere, i.e. the standard basis of R n . 

Then, the Weyl equation (J2]) takes the form 

s + { £ ft (<&. + «&*) « (4? + 4?) - 1 £ ft («S u - 410 ° (4? - 4? 

fc=i fc=i 

+{ £ °* 4l } ® (4? - 4?) - * I ® (el? - eg) - x I ® (eg + eg) - <y I ® (eg - eg 
fe=i 

n 

E^e[ n) ®(5ieg + <72eg)=0. (20) 
To proceed with our computation we exploit the following properties 

4°4?=^4 n \ efe^=5 jk t ] . (21) 
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With the use of the latter identities, and after setting 



v t (r, s) = h t (r, s) g^r, s), w t (r, s) = hi(r, s) g 2 (r, s), 



(22) 



equation (1201) is equivalent to the following first-order system of differential equations 



where k — 1, 2, . . . , n — 1. Here, Vi and Wi for i = 1, . . . , n are the total derivatives of the 
functions f«(r, s) and Wi(r, s) with respect to the argument s. 

Solving the aforementioned system is the first step in reconstructing the solution of the 
self-dual Yang-Mills equation from Nahm data. Then the problem of recovering the Higgs 
field and gauge potentials given by ([6]) is linear. However, solving the system (123!) is a rather 
technically difficult task. In the next section, we show how the existing technical difficulties 
can be overcome by connecting the corresponding full Weyl equation (T5]) with the isotropic 
Heisenberg model. 



The Weyl equation (T5]) can be identified as a Hamiltonian system. In particular, one can ob- 
serve that ([2]) can be described by a Hamiltonian containing some bulk spin-spin interaction 
and a boundary term. 

For simplicity, let us focus on the SU(3) case. The corresponding results can be extended 
in all other cases; however, not in a straightforward manner. After implementing the simple 
Nahm data of the minimal symmetry breaking fl3]), the corresponding full Weyl equation ([2]) 
can be written in the form 




v k +i - -C k w k - ( — a k+ i + zj v k+1 - (x - iy) w k+1 = 0, 
w k C k v k+1 + ( a k + z ) w k - (x + iy) v k = 0, 



(23) 



3 Weyl Equation & Open Spin Chains 




(24) 
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where 



H = — (cr x ® a x + a y ® a y + a z ® cr^) + x I ® a x + y I <g> a y + z I <g> cr^ (25) 

^5 



which is nothing else but the open XXX (Heisenberg) spin chain. The first term in (1231) 
corresponds to the spin-spin interaction; while the rest descibes the boundary interaction. 
This is a quantum integrable model and the corresponding Hamiltonian is immediately 
obtained from the so-called open transfer matrix. We shall not provide the details of this 
construction in the present investigation, but we refer the interested reader to the original 
work discussed in [T3] . 

The boundary part of the XXX Hamiltonian (1231) can be explicitly expressed as 

/ 1 tan0e _i * \ 
xo x + y<Jy + zo z = z . (26) 

V tanfle** -1 I 

It is clear that there is an one to one correspondence between the coordinates and the 
boundary parameters of the open XXX model. 

The key point in solving the generic problem described by ( 1241) and ( |25l) is the following 
observation: A similarity transformation U (see, e.g. [IT]) exist, which diagonalizes the 
boundary contribution in (1231) . but leaves the bulk spin-spin interaction invariant. More 
precisely, let 

/ 1 + cos 9 sin 9 e~^ \ 

U=[ . (27) 

\ (l-cos0)e** -sin^ I 

Then it is straightforward to show that, the boundary term is transformed as 

1 tan 9 . , . , _ 

u\ ., \u- l = —A j» 




tan0 -1 
while the bulk spin-spin interaction term remains unaffected: 

U ® U (cr x ® o x + o y ® a y + o z ® o z ) U^ 1 ® U^ 1 = cr x ® a x + a y ® o y + a z ® a z . (29) 

When the transformation ( 1271) acts to the full Weyl equation ( 1241) and ( |25l) . one gets: 

£/ ® C/ ( 4- ~ H J v = =^ I -T- - ;r- ^ ® °"i - r 1 ® <r z J ([/ ® E7) v = (30) 
\ds J \ds 2s J 
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which is nothing else but the Weyl equation with only diagonal boundary terms. In p.], so- 
lutions of the aforementioned equation have been obtained. Therefore, the inverse similarity 
transformation can be performed in order to solve explicitly the full non- diagonal problem. 

Let v to be the solution of the diagonal problem [1]. Then the solutions v of the full 
problem can be obtained from v due to the relation 

v = (U- 1 g> U- 1 ) v . (31) 

Let us also briefly consider the simplest case, i.e. the SU{2) BPS monopole. In this 
case, the corresponding Hamiltonian is just the boundary term ( 1261) . This Hamiltonian 
emerges from the open transfer matrix as an 1-site open XXX chain. After implementing 
the similarity transformation (127)1 one gets: 

u(^-n)v = ^ (-^-raAUv^O. (32) 

Therefore, the vector v in terms of the solution of the diagonal problem is of the form 

v = f/- 1 v . (33) 

In general, for the SU(n + 1) case the corresponding Hamiltonian has the form: 

U = ^-Ti ® <Ti + o-j. (34) 

Assume that there exist a transformation U (ie. nxn matrix) such that the bulk interaction 
term is 

U®U(Ti® o-i) W l <8> U~ l = (r; g> (Ji) . (35) 

Then, the aforementioned procedure can be applied in order to obtain all the solutions of 
the general SU(n + 1) case. The generic solution of the Weyl equation will then be 

v= {U- l ®lJ- 1 ) v , (36) 

where v the solution of the diagonal case found in pp. Hence, the main question to be 
investigated is the derivation of the transformation U. Extending our work on the 577(2) 
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and S77(3) it is obvious to see that the n x n matrix U in the general case SU(n + 1) is of 
the form 

u- 1 = (/, (i + If D^^. • • ■ , U + 1?! 2 )^^) < 37 > 

where / is an n-component column vector which is a holomoprhic function of £ = tan(#/2)e ?<?i 
of degree (n — 1), the so-called harmonic maps from S 2 to CP n . These spherically symmetric 
maps into CP n are given by 




/ (fo fj /„ ,)'• where /, . (38) 



and ( n ~ ) denote the binomial coefficients. It can be shown that these maps, which describe 
the instanton and non-instanton solutions of the CP n model [13] . are spherically symmetric 
in the sense that a rotation in R 3 exists which is realized as Mobius transformation. 

For n = 2 these are all the finite action solutions, but for n > 2 there are other non- 
instanton solutions. These can be described by introducing the operator P + defined by its 
action on any vector / 6 C" as 

P+f = drf - fif ^ f) (39) 

and then define further vectors P£/ by induction: P£/ = P + (P^" 1 f) . 

To proceed further we note the following useful properties of Pj?f when / is holomorphic: 



(Plf)ip l + f = 0, k^l 



I pk f 12 / pk-l f \ pk f 

a- (plf) = -p k - l f 'vf , a ( + J ) = V 7 ■ (40) 



These properties either follow directly from the definition of P + or are easy to prove [15J. 
Applying P + a total of n — 1 times to a holomorphic vector gives an anti-holomorphic vector, 
so that a further application of P + gives the zero vector. The modulus of the corresponding 
vector P?_f for / of the above form is 

\Plf\ 2 = a (i + l^l 2 )"- 2 ^ 1 , k > 1 (41) 

where a depends on n and k [16] ; while from construction \ f\ 2 = (1 + |£| 2 ) n_1 . Remark: In 
order our matrix U to leave the bulk interaction unaffected, ie. to satisfy the equation (I33|) . 
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the norm of each harmonic map should be equal to (1 + |£| 2 ) n_1 . Thus we have to normalize 
analogously and obtain ( 137|) . (See also, Appendix, for some specific examples). 

Finally, one needs to implement the orthonormality conditions to the solutions v(x, s). 
The idea is to pick n + 1 distinct solutions and require that they satisfy the orthonormality 
relations: 

<v i ,v i >= 5 i:j . (42) 

This way one may eventually identify the Higgs field and gauge potentials for the generic 
situation via In the next section two specific examples, displaying explicitly how the 
proposed methodology works, are discussed in detail. 



4 Explicit Examples 



4.1 577(2) Case 

Let us first consider the simplest case within the frame we described above, i.e. the SU{2) 
gauge. Recall that, the solution of the diagonal Weyl equation is of the form 



v 



92(r,s) 



(43) 



a(r) e rs 
b(r) e~ rt 

Then, the corresponding solution of the full Weyl equation, that is after performing the 
inverse similarity transformation ( 1331) . is given by 

I / - [a(r) e rs + b(r) e~ rs e~^] sin 9 

a(r) e rs (cos# - 1) e<* + b(r) e~ rs (cosfl + 1) 

with constants a,- and 



2sin# 

Next, chose two vectors Vj for % = 1, 2 to be of the form 
such that the set {vi, V2} defines an orthonormal system. One may readily extract that 



(44) 



< V i; Vj > 



1 f 1 

-^5— / \at a 7 - (1 - cos 9) e 2rs + 6< b j (cos 9 + 1) e~ 2rs ] ds 
m z 9 J_i J 



2 sin 

sinh(2r 
2r~ 



Cb^ (Xi j 



+ 



bjbj 
'1-u) 



(45) 
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using the variables fl9]). Then, from the orthonormality conditions (1421) the form of the 
constants cij and 6, which are now functions of the variables (r, u, ip) (up to a gauge choice) 
can be obtained. In fact, ( PHI) becomes ( IT5l) by setting 



«2 



sinh 2r 



61 = - 



sinh 2r 



(46) 



sinh2r ~ V sinh2r y ^ 

Finally, the associated Higgs field and gauge potentials are given by ( TTB"j) . For example, the 
elements of the Higgs field are of the form 

• 1 



a i a j e 2rs _|_ jjj bj g-2rs 



2J_ X L(l + «) (!-«) 
if — sinh 2r + 2r cosh 2r 



'/> 



2r 2 



(47) 



4.2 SU{3) Case 



The next non-trivial case, which can be easily treated is the SU(3) gauge. Recall that the 
solution of the diagonal case has the form 



' Vi(r, s) > 



v 



hi(r,s) 
h 2 (r, s) 



9i(r, s) 
Skir, s) 



(48) 



wi(r,s) 
v 2 (r, s) 
\ w 2 (r,s) J 

where Vi = gi hi and Wi — g 2 hi for i — 1, 2. The entries Ui and Wi are known explicitly due 
toP: 



v\ = a(r) \fs e r 



V2 = W+ ( r- — I 11. 



wi = — , 
s 

w 2 = b(r) \fse~ 



(49) 



where 



W = c(r)M[ -,l;2rs ] = c(r) 



sinh(rs) — rse 



(50) 
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M(k, m; z) is the Whittaker function of first type (see also pQ) and W denotes the derivative 
of the function W with respect to s. 

After performing the inverse similarity transformation the solution of the full Weyl equa- 
tion becomes 



1 



{hi + h 2 e-**) sin 9 



4sin 9\ hi(cos9-l) e i * + h 2 {cos6 + l) I \ gi (cos 9 - 1) + g % (cos 9 + 1) 



(51) 



As before, we choose three distinct vectors Vj for j — 1, 2, 3 of the form ( I5TT) . Then, the 
associated solutions v^\w^ are given in terms of the constants a,j,bj and Cj by ()49l) and 
( 150]) . Again, one may readily extract that their norm is equal to 



4sin 2 # 



(COS^-I) 2 M (j ) (COS^ + I) 2 



sin 2 6 



sin 2 6 



*W (i) i *(0 0) i *(*) (i) 



'/- 



1 

16 



a,- a 



9o 



'"V (l + uf 



^ 2 + CiCj ^ 2 



(52) 



r 2 (1 -uf " J 9r\4>\ 
where qo = 1 + (—1 + 6r) e _6r , gi = 1 — (1 + 6r) e 6r and (r, w, ■?/>) are the parameters defined 
in ([9]). Note that, for convenience the limits of integration have been shifted in this case to 
and 3. 

Requiring the orthonormality conditions (j4"2"|) one may fix, not uniquely, the constants 
which are now functions of (r,u,ip) and determine the vectors Vj. It is straightforward to 
show that the constants aj, bj, and Cj are of the form: 



4r 



a i 

a>2 



4r 



{u + 1) e* 
{u + 1) e* 



V 3( ?o 



4r 



v3gi 

4r 



4r 



(1 - u) e*, ci : 
(l-n)e^ + ¥, C2 , 

(l-u)e^ + f, c 3 



4y/r 



a/ 3 <?o q~i 
V 3 <?o <?i 



(53) 



where a, /3, and 7 are real free parameters, which for simplicity may be set equal to zero. 

Then, the Higgs field and gauge potentials can be recovered from ([6]). In particular, the 
elements of the Higgs field are equal to 



4|^l 2 



(s-2) 



{u-lf (i) (n + 1) 2 . *(i)Jj) . „„*(*)„„(j) 



u l «l + 



ds 
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2^ e 6r 



k bj 



l-2r+ (l-4r + 6r 2 ) 

4^ y ~ (1 -u) 2 

-36r 2 - 4 + (2 - 3r) e 6r + (2 + 3r) e~ 6r \ 
18r J 



-1 + 2r + (l +4r + 6r 2 ) e -6 ^ 



(54) 



5 Conclusions 

SU(n + 1) spherically symmetric monopole solutions of the full Weyl equation in the case 
of minimal symmetry breaking were obtained. This was done by implementing azimuthal 
dependence to the solutions found in [1] via a suitable similarity transformation. More pre- 
cisely, the existence of a spin chain-like (XXX chain) part with a bulk spin-spin interaction 
and a boundary term within the Weyl equation was exploited. In this context, a similar- 
ity transformation exists that turns the diagonal boundary terms to generic non-diagonal 
ones but, leaves the bulk interaction invariant [17]. This is precisely the transformation one 
utilizes to implement the azimuthal dependence to the spherical symmetric solutions. 

In this paper, a particular case that involves the simplest Nahm data is considered and 
thus, the isotropic Heisenberg model (XXX chain) is involved. For generic Nahm data the 
situation becomes more complicated, and the anisotropic Heisenerg model (XX Z or XYZ 
chain) is involved. In the general situation, such a boundary similarity transformation that 
leaves the bulk part unaffected does not exist; thus, more sophisticated methods need to be 
employed (see, for example, References [HI HI]). This is an intriguing issue, that we hope 
to address in future investigations. 



A Appendix 

We give two explicit examples of the transformation U given by (l37j) for the SU(4) and 
SU(5) case. 
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SU (4) Case: Here n = 3 and our formulation leads to 



/ 1 f X 



ZT 1 



\f V2( 1 
SU(5) Case: Here n = 4 and our formulation leads to 



u- 1 



vse — e ciei 2 — 2) -(2iei 2 -i) -vse 



/ 



(55) 



(56) 
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